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^ ■ ABSTRACT 

The Lyman decrement associated with the cumulative effect of H i in QSO absorption systems along 
,; the line of sight provides a distinctive feature for identifying galaxies at z > 2.5. Color criteria, which are 
J — ' sensitive to the presence of a Lyman-continuum break superposed on an otherwise flat UV spectrum, have 
T-H ' been shown, through Keck spectroscopy, to successfully identify a substantial population of star-forming 
^ . galaxies at 3 < z < 3.5 (Steidel et al. 1996a). Such objects have proven surprisingly elusive in field- 
\^ ' galaxy redshift surveys; quantifying their surface density and morphology is crucial for determining how 
Q\ ■ and when galaxies formed. The Hubble Deep Field (HDF) observations offer the opportunity to exploit 
[ the ubiquitous effect of intergalactic absorption and obtain useful statistical constraints on the redshift 
Qh' distribution of galaxies considerably fainter than current spectroscopic limits. We model the H i cosmic 
Q . opacity as a function of redshift, including scattering in resonant lines of the Lyman series and Lyman- 
' continuum absorption, and use stellar population synthesis models with a wide variety of ages, metallicities, 
^ dust contents, and redshifts, to derive color selection criteria that provide a robust separation between high 
• • , redshift and low redshift galaxies. From the HDF images we construct a sample of star-forming galaxies at 
' 2 < z < 4.5. While none of the ~ 60 objects in the HDF having known Keck/LRIS spectroscopic redshifts 
^ ' in the range < z < 1.4 is found to contaminate our high-redshift sample, our color criteria are able to 
^ ', efficiently select the 2.6 ^ z < 3.2 galaxies identified by Steidel et al. (1996b). 

The ultraviolet (and blue) dropout technique opens up the possibility of investigating cosmic star and 
element formation in the early universe. We set a lower-limit to the ejection rate of heavy elements per unit 
comoving volume from Type II supernovae at {z) = 2.75 of ~ 3.6 x 10~^ Mq yr~^ Mpc~^ (for qq = 0.5 and 
Hq = 50kms~^ Mpc~^), which is 3 times higher than the local value, but still 4 times lower than the rate 
observed at z ~ 1. At (z) = 4, our lower limit to the cosmic metal ejection rate is ~ 3 times lower than the 
(z) = 2.75 value. We discuss the implications of these results on models of galaxy formation, and on the 
chemical enrichment and ionization history of the intergalactic medium. 

Subject headings: cosmology: observations - galaxies: evolution - intergalactic medium - quasars: absorption 
lines - ultraviolet: galaxies 
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1. INTRODUCTION 

Much observing time has been devoted in the past few years to the problem of the detection of galaxies 
at high redshift, as it is anticipated that any knowledge of their early luminosity and color evolution will 
put constraints on the history of structure and metal formation. While it has become clear that blank-sky 
surveys for strong Lya -emitting primeval galaxies are not particularly efficient (e.g., Pritchet & Hartwick 
1990; Djorgovski &; Thompson 1992), the method of obtaining multicolor broadband observations of the 
emitter's rest-frame UV stellar continuum has been successfully applied to detect galaxies at cosmological 
distances (Steidel et al. 1996a). 

Ground-based observations have used color techniques which are sensitive to the presence of a Lyman- 
continuum break to identify or set limits on the number of high redshift galaxies (Stcidcl & Hamilton 1992, 
1993). The signature of a distant, star-forming galaxy in such surveys is a very red U—B color, combined with 
colors in longer-wavelength filters that are much bluer. The early work by Guhathakurta et al. (1990) showed 
that Lyman-break objects do not dominate the galaxy counts at faint apparent magnitudes. Assuming that 
young, star-forming galaxies have flat spectra longward of the Lyman-break, Steidel &; Hamilton (1992, 1993), 
and Steidel, Pettini, & Hamilton (1995) have designed and used a custom set of broad-band filters to select 
high-z galaxies in the fields of distant QSOs. Recent deep spectroscopy with the W. M. Keck telescope 
has shown the high efficiency of such color-selection technique: Steidel et al. (1996a) have identified 17 
Lyman-break star-forming galaxies with redshifts 3.0 ^ z ^ 3.5, which are most likely the progenitors of the 
present-day bright spirals and ellipticals observed during the assembly of their cores (Giavalisco, Steidel, & 
Macchetto 1996). 

With its high spatial resolution, the Hubble Space Telescope (HST) offers the opportunity to study the 
faint galaxy population in unprecedented detail. In particular, the images of the Hubble Deep Field (HDF), 
obtained with the Wide Field Planetary Camera (WFPC-2) during December 1995, represent the deepest 
optical imaging survey undertaken so far (Williams et al. 1996), reaching S-u limiting AB magnitudes of 
roughly 27.7, 28.6, 29.0, and 28.4 (for an aperture area of 0.2 square arcsec) in the F300W, F450W, F606W, 
and F814W bandpasses (the number corresponds to the central wavelength in nm), respectively. In this 
paper, we describe selection criteria for this filter system that can be used to identify a sample of likely 
precursors to present day galaxies at 2; > 2, with little contamination from low redshift galaxies. Our 
strategy is designed to exploit the effect of the increasing opacity of the intergalactic medium (IGM) at high 
redshifts, as indicated by the plethora of absorption lines seen in the spectra of background quasars: enough 
neutral hydrogen is known to be contained in the clumps of highly ionized gas which form the Lya forest, 
and in the metal-line absorption systems associated with the halo regions of intervening bright galaxies, to 
significantly attenuate the UV flux from distant sources. We show that the cumulative effect of H i in QSO 
absorption systems along the line of sight provides a distinctive feature for identifying star-forming galaxies 
in the HDF. Although other spectral features, such as the 4000 and 912 A breaks which characterize the 
integrated spectra of stellar populations, with the latter possibly enhanced by self-absorption from interstellar 
gas within the galaxy itself, may also be used as tracers of redshifts, the model predictions of their magnitude 
are sensitive to the unknown physical and evolutionary state of the galaxy, i.e., its star formation history, 
age, and H i distribution, and hence are subject to substantial uncertainties. By contrast, the "reddening" 
effect due to atomic processes in cosmological distributed QSO absorption systems is ubiquitous, quite strong 
for z > 2.5, and can be reliably taken into account. Although stochastic in nature, r.m.s. fluctuations away 
from the mean opacity are bound to be modest in most situations, due to the broadband nature of the 
adopted filter set. 

The technique to analyze the HDF images we develop in this paper is based on the theoretical study 
by Madau (1995) (see also Yoshii &: Peterson 1994), and is an extension of the Lyman-break color criterion 
developed by Steidel & Hamilton (1992, 1993), and Steidel, Pettini, & Hamilton (1995). We extend this 
technique to the HDF filter system, and perform several experiments with simulated and real galaxy spectra 
to define effective selection criteria for high redshift galaxies, and identify possible sources of contamination 
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from low redshifts. By computing colors for an extremely wide range of model galaxy spectra, we are able 
to tune the criteria to provide what we believe are largely uncontaminated samples of star-forming galaxies 
in the redshift ranges 2 < z < 3.5 and 3.5 < z < 4.5, and use these samples to estimate the integrated star 
formation and metal ejection rates at those redshifts. We show how follow-up spectroscopy with Keck/LRIS 
by Steidcl et al. (1996b), the Hawaii+Caltcch group (Cohen et al. 1996), and the Berkeley HDF group 
(Moustakas, Zepf, & Davis 1996) supports the efficiency of our color selection techniques. Future work 
on high redshift galaxies in the HDF will investigate the luminosity function, study the morphologies, and 
compare the numbers and properties of high-z candidates to the predictions of a variety of galaxy-evolution 
models. 

Throughout this paper, unless otherwise stated we shall adopt a flat cosmology with go = 0-5 and 
Ho = 50kms"^ Mpc"^. 

2. INTERGALACTIC ATTENUATION 

In this section we will briefly review the theory of the propagation of UV radiation through a clumpy 
universe, following Madau (1995). 

2.1. Basic Equations 

Let -L(^'em) be the specific power emitted with frequency fem by a source at redshift Zem- The mean 
specific flux observed at Earth is 

{fM)=^-^^^^^^{e-l, (1) 

L 

where fobs = t'em/(l + -^em), t^L IS the luminosity distance to z^m, and the average transmission over all lines 
of sight is, assuming Poisson-distributed clouds. 

Here, is the optical depth through an individual cloud at frequency u = fobs(l + -2), and {d'^N /dN-mdz) 
is the redshift and column density distribution of absorbers along the path. An "effective" optical depth of 
a clumpy medium can be defined as Teff- = — ln((e~'^)). 

Along with resonant line scattering from Lyct, /?, 7, and higher order members, we include photoelectric 
absorption from H i in the Lya forest clouds and Lyman-limit systems along the line of sight. Since the 
bluest filter used for the HDF observations is centered at 3000 A and has FWHM~ 800 A, galaxies will only 
be subject to H i cosmological attenuation from material located at z > (2600/1216) — 1 = 1.1. 

2.2. Line Blanketing & Continuum Absorption 

For A/3(l -I- Zem) < Aobs < '^a(l + -^em), where Afj = 1216 A and = 1026 A, a galaxy's continuum 
intensity is attenuated by the combined blanketing effect of many Lya forest absorption lines, with effective 
opacity 

/ , \ 3-46 

Teff = 0.0036 f^j , (3) 

(Press, Rybicki, k, Schneider 1993). Hence, line blanketing from Lya alone will produce > Imag of attenu- 
ation in the continuum spectrum shortward of 6000 A of a galaxy at ^ 4. 

When Aobs < -^/?(1 + -^em), a significant contribution to the blanketing opacity comes from the higher 
order lines of the Lyman series. A standard curve of growth analysis has been applied to numerically 
compute the attenuation expected from line blanketing of Ly/3, 7, 5 plus 13 higher order members. In the 
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wavelength range Ai_|_i(l + Zem) < Aobs < Aj(l + Zem), the total optical depth can be written as the sum of 
the contributions from the j ^ 1 transitions, 



where Aj = (1.7 x lO'^, 1.2 x 10-3,9.3 x lO""*), and A^- = (1026,973, 950 A) for Ly/3, 7, and 6, respectively. 
It is worth noting that heavy element absorbers make a negligible contribution to the blanketing optical 
depth at high redshifts, as this is dominated by those lines which lie at the transition between the linear 
and the flat part of the curve of growth, i.e., with A^hi ~ 10^3-6 cm"^ in the case of Lya . 

Continuum absorption by H i affects photons observed at Aobs < Al(1 + Zem), where Al = 912 A is the 
Lyman limit. An approximate (within 5%) integration of equation (2) yields for the effective photoelectric 
optical depth along the line of sight: 



r^ff = 0.25xl{xl:^' - xl-^') + 9Axl-Hx'^i' - - 0.7xUx-'-'' - x-^^^) - 0m3{xl^' - xl'^^ (5) 



where Xc = (Aobs/AL) for Aobs > ^l, and x^m = 1 + z^m- The first term on the right-hand side represents 
the approximate contribution of Lyman-a clouds, the others are due to Lyman-limit systems. Absorbers 
with iVni ~ 10^^ cm~^ dominate the cosmic continuum opacity.^ 

While our formula for the line blanketing optical depth is entirely consistent with the flux deficits 
Da and Db observed in the spectra of QSOs below their Lya and Ly/3 emission features (e.g., Schneider, 
Schmidt, & Gunn 1991), the derived continuum opacity is subject to significant uncertainties, as very limited 
information exist on optically thin absorbers with 10^^ < A^hi ^ 10^^ cm~^. 



Figure la shows the characteristic staircase profile (cf M0ller k, Jakobsen 1990) of the mean cosmic 
transmission (e""^) for a source at Zem = 3, 4, and 5, as a function of observed wavelength. Also plotted are 
the response functions of the photometric system we shall adopt in our discussion. This consists of the four 
broad ifST passbands F300W, F450W, F606W, and F814W (roughly UBVI). These filters efficiently cover 
most of the bandpass accessible to the HST WFPC-2. 

It is clear that the observed broadband colors of cosmological distant objects will be strongly "reddened" 
by continuum absorption and blanketing of discrete absorption lines which move into and out of the color 
passbands with changing z. To be quantitative, we must take into account that the mean transmission 
observed is not (e""^) but rather an average over the bandpass. 



where r(A) is the normalized transmittance of the relevant filter. In the following, we shall measure the mean 

integrated intergalactic attenuation in magnitudes. Am = — 1.086 In Q. Figure Ih shows A[/3oo, AB450, 
^^606) and A/814, the observed magnitude increases at the corresponding bandpass due to intervening 
absorption, as a function of emission redshift in the redshift range 1.5 < z < 5. [These increments must be 
added to the term — 1.086 In (1 + Zem) to get the standard K-correction for a flat emitted spectrum.] 

We are now able to assess how intergalactic absorption will modify the intrinsic photometric properties 
of galaxies at high redshift. At Zem ~ 1-5, Lya line blanketing starts to cause a small apparent depression 
in the C/300-band continuum. As higher order lines of the Lyman series move into the bandpass, At/soo 

^ The He i contribution to the attenuation is negligible in the case of a QSO-dominated ionizing back- 
ground, while He 11 absorption on the way must be included if (1 + Zem) > Aobs/228A (e.g., Madau 1992). 




(4) 



2.3. Cosmic Transmission 




(6) 
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iiicrcases rapidly until, at Zem = 2.3, sources in the background appear, on average, ~ 0.7 mag fainter in 
the ultraviolet. By = 3, Lyman-continuum absorption from forest clouds and Lyman-limit systems is so 
severe that AUsqq « 2.8 mag (of which « 1.9 mag are due to Lya clouds alone): galaxies will either appear 
very red in L^soo — -B450 or will drop out of the F300W image altogether.^ Similarly, the Lya line progressively 
enters the F450W band beyond ^ 2.7. At Zem = 3.8, line blanketing from the Lyman scries produces 
A5450 ~ 1.2 mag. By = 4.4, Lyman-continuum absorption contributes significantly to the total opacity, 
and Ai?45o exceeds 3.3 mag. The F606W bandpass is only weakly affected by intergalactic absorption for 
Zem ^ 4. Hence, again, star-forming galaxies at Zem ~ 4 will either appear very red in B450 — Veoe or will 
effectively be undetectable in the blue band. By = 4.6, line blanketing from the Lyman series produces 
AVgoe ~ Iniag. As H i absorption does not affect the F814W passband for Zem ^ 5.7, A/814 = in the 
redshift range of interest here. 

Although Figure lb refers to the average accumulated absorption only, r.m.s. fluctuations away from 
the mean transmission are predicted to be relatively small after integration over broad bandpasses (Press 
et al. 1993). For example, the observed i?45o-band continuum of a galaxy at Zem ^ 3 is attenuated by 
the blanketing of more than 90 Lya lines with W > 0.3A along the light travel path. At these redshifts, 
excursions away from the average curve will be the largest in the F300W filter, due to the contribution to 
the photoelectric opacity by the rarer, optically thick Lyman-limit systems. Even in this case, however, with 
dN/dz « 0.27(1 -|- zY'^^ absorbers per unit redshift with A^hi > 1-6 x lO^^cm"^ (Storrie-Lombardi et al. 
1994), as many as AN ~ 1.5 Lyman-limit systems along a random path to Zem ^ 2.7 would be expected, on 
average, to photoelectrically absorb photons in the ultraviolet bandpass. As the probability of intersecting 
at least one metal system is large, 1 — e~^^ ~ 80%, we expect only one out of five star-forming galaxies at 
z > 2.7 to be effectively detectable in the ultraviolet bandpass. 



One of the major objectives of the HDF observations is to determine the morphology, spectral energy 
distribution, and number density of galaxies actively forming stars at high-z, the nature of which is crucial to 
our understanding of galaxy evolution. In this section we shall show how we can use intergalactic absorption 
to identify high redshift galaxies in broadband multicolor surveys. Colors can be more readily interpreted 
using the AB-magnitude system (Oke 1974) 



where f{u)e~'^'ff is the mean incident power measured in ergs cm~^ s~^ Hz~^, and T{v) is the normal- 
ized transmittance of the relevant filter. We put our magnitudes directly into the AB system, so that 
[/300 = -B450 = Veoe = -^814 corresponds to a spectrum for which f{u)e~'^''ff is constant. The observed 
colors of distant objects are determined by the intrinsic spectral energy distributions (SEDs) of their stellar 
populations, absorption and scattering of photons in their interstellar media, and stochastic absorption due 
to intervening H i clouds. Because we do not know when galaxies started forming stars, and over what 
timescale they formed them, it is impossible to predict their energy distribution in detail. Furthermore, it is 
not necessarily the case that the faintest galaxies are the most distant. Thus our analysis must not rely on 
fluxes or apparent magnitudes to give redshift information. It is, however, possible to outline in a general 
way the range of plausible colors for star- forming galaxies at high redshift, and to identify combinations of 
age and dust reddening in a stellar population at lower redshift that may conspire to mimic a high-z galaxy. 



^ Note that AU300 flattens for Zem ^ 3.6 because of a significant "red-leak" of the F300W filter above 
6700 A. 



3. THE COLORS OF HIGH REDSHIFT GALAXIES 




(7) 



3.1. The Spectral Energy Distribution of Star-Forming Galaxies 
For illustrative purposes, we analyze here two simple examples in detail: 
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(1) We compute the ultraviolet spectrum of a star-forming galaxy using the isochrone synthesis spectral 
evolutionary code of Bruzual &: Chariot (1993). The model has solar metallicity, a constant star formation 
rate at age 0.3 Gyr, a Salpeter (1955) initial mass function (IMF), (p{m) oc m~^ '^^, with lower and upper 
cutoffs of 0.1 and 125 M©, and completely ignores the effects of H i in the local interstellar medium on the 
transfer of ionizing photons. Two intrinsic discontinuities, about 1.4 mag across the rest-frame Lyman limit 
and 0.5 mag across the Balmer decrement, are characteristic of the most prominent stars; 

(2) We take as fiducial the spectra of the four brightest starburst galaxies observed by the Hopkins 
Ultraviolet Telescope (HUT) on Astro-2: NGC3310, 4214, 5236 and 5253 (Ferguson et al. 1996a). The HUT 
spectra cover the wavelength range 912 < A < 1820 A (Kruk et al. 1995). The cutoff at short wavelengths 
is due to galactic H i . For wavelengths longward of 1820 A, we have patched on the lUE spectra from the 
Kinney et al. (1993) atlas, renormalized to the HUT flux in the range 1350 — 1700 A. The renormalization 
is required because of the different sizes of the HUT and lUE apertures. Note that the starburst galaxy 
intrinsic colors are significantly redder than the ones derived from the population synthesis model because of 
the presence of a large Lyman discontinuity and of dust-reddening. i^C/T spectra of somewhat more distant 
starbursts galaxies suggest that the intrinsic Lyman break is typically at least a factor of 5 (Leitherer et 
al. 1995). Analysis of the NGC5236 HUT SED suggests that a young (~ 5 x 10^ year old) population is 
responsible for the emission, with an intervening extinction E{B — V) ~ 0.2 necessary to fit the detailed SED 
(Ferguson et al. 1996a). Such a modest extinction appears fairly typical of starburst galaxies (even those of 
high metallicity) . In a comprehensive study of 39 starburst galaxies, Calzetti, Kinney, & Storchi-Bergmann 
(1994) found values of E{B — V) ranging from to 0.8. 

Figure 2a displays the predicted U^qq — i?45o, ^450 ~ ^06 1 ^^nd ^06 " -^814 colors versus emission redshift 
of our synthetic IGM-attenuated galaxy, together with the unattenuated values. The intrinsic spectral shape 
has been kept constant with cosmic time. As expected, the observed colors redden quite sharply with redshift 
due to intergalactic absorption. For example, from = 2.5 to Zem = 3, the observed [/300 — -B450 color of 
the Bruzual spectrum increases from 1.8 to 3.6 mag. In the absence of any absorption break, the last value 
would correspond to a spectral energy distribution /(v) oc with a = 87 The observed flux decrements 
between 4500 and 6060 A and between 6060 and 8140 A similarly probe the universe at higher and higher 
redshifts, as i?45o — Vqqq increases from 0.8 at Zem = 3.5 to 1.7 mag at Zem = 4, while Vgoe ~ -^814 increases 
from 1.0 at Zem = 4.5 to 1.9 mag at Zem = 5. A similar plot is shown in Figure 2b for the observed HUT-IUE 
starburst galaxy spectra. The main point to note here is how intrinsically blue and red objects follow similar 
paths in the color-redshift plane, and that, even in the presence of dust-reddening, the large magnitude 
jumps between adjacent bandpasses are largely caused by the known sources of intergalactic H i opacity - the 
Lya clouds and Lyman-limit systems - rather than by the assumed galaxy SED. 

3.2. Color Selection Criteria 

The problem of identifying high redshift galaxies is complicated by the fact that the present-day colors 
and spectra of galaxies allow for a very wide variety of star-forming histories. This evolution is coupled with 
chemical evolution of the stars and the ISM and changes in the amount of dust and its distribution relative 
to the stars. However, galaxy spectra are not entirely arbitrary. Flux is usually a slowly varying function of 
wavelength, with a few spectral discontinuities (e.g., the 4000 and 912 A breaks) at specific locations. Thus 
it should be possible to construct fairly robust selection criteria that will exploit the combined effect of the 
intrinsic Lyman edge in galaxies and the opacity of intergalactic neutral hydrogen to separate high redshift 
from low redshift objects. 

In this section, we shall define criteria appropriate to the HDF bandpasses. To derive a robust color 
selection technique, we have computed HDF colors for 1612 synthetic spectra of galaxies, representing 



^ Above Zem ~ 3.5, J/aoo — -B450 becomes bluer, again an effect of the red-leak of the F300W filter. 
However, by then, intergalactic attenuation is so strong that high-z galaxies are effectively unobservable in 
the ultraviolet. 
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different ages, star formation histories, metallicities, and dust opacities. The stellar-evolutionary input to 
the models has been described in detail by Babul & Ferguson (1996). Briefly, synthetic spectra are computed 
using isochrone synthesis from the isochrones of Bertelli et al. (1994), and the model atmospheres of Kurucz 
(1992) and Clegg &: Middlemass (1987). The synthetic spectra agree well with those of Bressan, Chiosi, 
& Fagotto (1994), which were constructed using the same isochrones, and slightly less well with those of 
Bruzual & Chariot (1993), primarily due to their use of different isochrones. The uncertainties in population 
synthesis have been outlined in some detail by Chariot, Worthey, &; Bressan (1996). These uncertainties 
amount to differences in color of several tenths for the same stellar population computed by different codes. 
While they can cause very large discrepancies in the derived ages, such small variations in color are almost 
completely negligible for our purposes. 

Attenuation by dust has been included using the parametrized extinction laws of Pei (1992) for the 
Galaxy, the LMC and the SMC. Note that the use of the relatively "gray" extinction curve for starburst 
galaxies found by Calzetti et al. (1994) would imply a higher UV continuum than predicted on the basis 
of the optical continuum emission by the application of the standard extinction laws, hence strengthen the 
efficiency of our selection criteria. We ignore the effects of H i in the local interstellar medium of the galaxies, 
which will add to the intergalactic attenuation to make high-2; objects even redder. There are no emission 
lines included in the simulations. Given the very wide bandpasses of the HDF filters, we do not expect 
emission lines to significantly alter galaxy broadband colors. 

The colors for the models have been computed at 73 different redshifts spanning the interval 0.001 < 
z < 7, taking into account the effect of intergalactic attenuation. In this exercise, galaxies arc required 
to be younger than the age of the universe (for qq = 0.5, Hq = 50kms~^ Mpc"-*^), but otherwise age has 
been decoupled from redshift. The model grid is intended to span the range of plausible colors of real 
galaxies, but not to sample this range in a way that is tied to cosmological models. For example, in the 
real universe we expect the old, low redshift portion of a color-color diagram to be more densily populated 
than in our simulations. The adopted grid of ages, star formation tiniescales, extinctions Ab, and redshifts 
is summarized in Table 1. It is important to emphasize that our goal is to construct a robust set of color 
selection criteria that are largely independent of models of galaxy formation or the assumed values of the 
cosmological parameters. In this respect our approach is quite distinct from efforts to assign photometric 
redshifts to each galaxy, as presented by Lanzetta, Yahil, & Fernandez-Soto (1996) and Gwyn &; Hartwick 
(1996), for example. 

3.2.1. F300W Dropouts 

We are now in position to examine the colors of arbitrary galaxies as a function of redshift. Figure 
3 shows i/soo — -B450 vs. S450 — /814 for all of the synthetic spectra. The large points show objects with 
2 < z < 3.5. Galaxies in this redshift range predominantly occupy the top left portion of the plot because of 
the attenuation by the IGM and intrinsic extinction. Galaxies at lower redshift can have similar 6^300 — ^450 
colors, but they are typically either old or dusty, and are therefore red in -B450 — Isu as well. To identify 
likely star-forming galaxies in our desired redshift range, we require that they have (1) C/300 — -B450 > 1.3, 
(2) [/300 — > -B450 — ^814 + 1-2, and (3) i?45o — I8i4 < 1-5. These criteria isolate objects that have 
relatively blue colors in the optical, but a sharp drop into the UV. 

Figure 4 gives an estimate of the eflficiency of this technique. The solid histogram shows the fraction 
of the total sample of models that meet the selection criteria as a function of redshift. Even with the wide 
range of parameters in our galaxy models, no galaxies with z < 1.5 or z > 4 are selected. However, many 
of the galaxies in the redshift range 2 < 2; < 4 are missed as well. The galaxies that are missed tend to 
be redder than the selection line in B450 — hiA, either because they are relatively old, or because they are 
highly reddened by dust. The dotted histogram in Figure 4 shows the selection efficiency for galaxies with 
ages restricted to be less than 10^ yr and extinctions Ab < 2. The selection criteria appear to be extremely 
efficient, recovering roughly 90% of the galaxies in the redshift range 2 < z < 3.5. Of course the true 
efficiency depends on how close our models are to real protogalaxies, and on how many old or very dusty 
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galaxies there are in the real universe in that redshift range. 

3.2.2. F450W Dropouts 

Wc can play an identical game with the F450W dropouts. Figure 5 shows the colors of the model galaxies 
in the ^450 — ^06 vs. V606~^8i4 plane, with galaxies in the redshift range 3.5 < z < 4.5 highlighted, and with 
our proposed selection criteria drawn. In this case, there are objects outside the redshift range of interest 
that can have rather blue Vgoe ~ -^814 colors and red -B450 — Veoe colors, and hence potentially contaminate 
the sample. These peculiar colors obtain only for those galaxies that lie in the redshift range 2 < z < 3.5 and 
have a fair amount of internal extinction. The blue Vgoe ~ -^814 color is due to the 2200 A absorption feature 
shifting into the F814W bandpass. Our selection criteria are tuned to avoid including many of these galaxies. 
Specifically, we require our candidates to have (1) S450 — Vgoe > 1-5, (2) B450 — Veoe > 1-7(V606 —-^814) +0-7, 
(3) S450 - "^606 < 3.5(^606 - /814) + 1.5, and (4) Veoe - ^8i4 < 1-5. 

Figure 6 shows the fraction of galaxies recovered as a function of redshift for such a selection. Once 
again, the criteria appear to be very efficient at selecting relatively unreddened star-forming galaxies, this 
time in the redshift range 3.5 < 2; < 4.5. Neverthless, there is some danger of low redshift contamination 
for the F450W dropouts. Again, most of the contaminants are either relatively old or reddened galaxies at 
slightly lower redshift. The dashed histogram in Figure 6 shows the distribution of redshifts for the sample 
when galaxies are restricted to have either ages greater than 10^ yr, or Ab > 2. Again, the extent of the 
contamination in the real universe depends on the unknown distribution of dust contents and ages at these 
redshifts. 

4. HIGH REDSHIFT GALAXIES IN THE HUBBLE DEEP FIELD 

Using the color selection criteria described above, we shall build in this section a sample of galaxies at 

high- 2; selected from the Version 2 catalog of the Hubble Deep Field, which is based on the second release of the 
reduced HDF images. Details of the data reduction, image combination, source detection, and photometry 
are given by Williams et al. (1996). 

Figures 7a and 8a show color-color plots of the catalog data with our selection region drawn in as 
dashed lines. The catalog is based on the three WF chips only, and covers an area of 4.65 arcmin^. The 
source detection was carried out on the exposure-weighted sum of the F606W and F814W images, and 
the threshold was set to 4(t. In all cases, isophotal magnitudes and colors were used for photometry, with 
identical apertures (defined from a sum of the F606W and F814W images) applied to the images in all 
bandpasses. It is important to note that the FOCAS isophotal magnitudes are systematic underestimates of 
the true brightness of the galaxies. The magnitude error is difficult to quantify precisely as it depends on the 
surface brightness profiles of the galaxies in the image. Experiments on simulated data suggest that total 
magnitudes are typically 0.5 mag brighter than FOCAS isophotal magnitudes (Ferguson et al. 1996b). We 
have ignored this offset in our analysis, but note that it could increase our inferred luminosity densities and 
metal production rates by roughly a factor of 1.6. The catalog was truncated at magnitude limits sufficiently 
bright that Lyman breaks of the expected amplitude could be measured reliably. For the U BI selection, we 
considered only objects with i?45o < 26.79 and Vf,o6 < 28.0; for the BVI selection, the limiting threshold 
was set to Vgoe < 27.67. Galaxies undetected in F300W (F450W) - as defined by having signal-to-noise ratio 
< 1 inside the isophotal aperture - were assigned a la lower limit to their U^oo — -B450 (-B450 — Veoe) colors. 
In the figures, only objects with /814 > 21 are shown. For the UBI selection of F300W dropout objects, 
two bright stars (with I814 < 21) enter into the color selection window. Steidel et al. (1996a) have shown 
that subdwarf stars may have similar colors to those of 2; ~ 3 galaxies, and indeed Steidel et al. (1996b) 
confirm that one of the two bright stars in the HDF within this color selection window is indeed a foreground 
subdwarf star. None of the fainter F300W dropouts appears to be stellar. 

The complex, sometimes multi-component morphologies of faint galaxies in deep HST images provide 
a complication for anyone cataloging faint objects in images such as the HDF. Automated object detection 
routines may break up objects into several components if their light distribution has multiple peaks, and it 
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is not always clear whether a multi-component object should be considered as a single "parent" entity or 
split into separate "daughters." We chose to apply our selection criteria to the complete catalog of both 
parents and daughters, and then inspected all potential high redshift candidates visually. In a few cases, we 
retained "parents" as single objects (e.g., for the brightest F300W dropout, the elongated object known as 
C4-06 in the notation of Steidcl ct al. 1996b), while for others we chose the "daughters" when the parent 
seemed to comprise distinctly separate objects. Although we shall consider here only the total, summed 
light from all candidate Lyman break objects, the details of this splitting can be of some importance as the 
daughters may fall in or out of our sample. 

Ultimately, of course, the ability of our color criteria to select largely uncontaminated samples of star- 
forming galaxies at high-z must be confirmed by deep spectroscopy. In Figures 7b and 86 we plot the 
location in the color-color diagrams of the more than 60 galaxies in the HDF which have known Keck/LRIS 
spectroscopic redshifts (Steidel et al. 1996b; Cohen et al. 1996; Moustakas et al. 1996). While none of the 
\ow-z objects is found to contaminate our high-z selection regions, our UBI color criteria are able to select 
4 of the 5 objects spectroscopically confirmed by Steidel et al. (1996b) to be galaxies in the redshift range 
2.591 < z < 3.226. (Here, photometry from the Version 2 catalogs has been used, whereas Steidel et al. used 
Version 1 data.) Two out of the six known z > 2 objects in the HDF are not selected by our criteria. One (a 
z = 2.268 galaxy from the Cohen et al. sample) only just misses being included - its U — B color is slightly 
bluer than allowed by our color criteria. The second missed source (C4-09 from Steidel et al. ) is a more 
puzzling case - this object is a small diamond-shaped configuration of four distinct sub-components, each 
of which "disappears" in the F300W bandpass. Steidel et al. considered it as a single object and measured 
a redshift of 3.23. We fail to identify it as a candidate because it appears as a formal detection {2.7cr) at 
F300W in the Version 2 catalog, whereas Steidel et al. assigned only an upper limit to its [/300 flux. This 
points to some uncertainty in our current understanding of the noise properties of the F300W images, but we 
do not expect this to significantly impact the results presented here. In the Version 2 catalog, the brightest 
of the four sub-components of this small "quad" is individually classified as a dropout above our magnitude 
limit and with the appropriate colors, and is included in the analysis presented here. 

4.1. F300W Dropouts: Limits to the Luminosity Density at 2 < z < 3.5 

We have identified 69 F300W dropouts in the HDF images which satisfy the criteria established for 
star-forming galaxies at 2 < z < 3.5. Assuming this redshift interval has been uniformly probed, we derive 
a comovmg galaxy number density of 4.2 x 10"^ Mpc'^ (0.92 x 10"^ Mpc"^ for go = 0.05) at (z) = 2.75, 

about 10 times larger than the comoving space density of bright L > L^ present-day galaxies (Loveday et al. 
1992). The observed Veoe magnitudes of our sample yield a specific (comoving) emissivity at 1620 A equal 
to 

P1620 !^ 1.6 X 10^6 ergs s"^ Hz"^ Mpc"^ (8) 

(9.0 X 10^^ ergss-i Hz"^ Mpc'^ for qo = 0.05). The faintest F300W dropout has a magnitude of Vooe = 26.7, 
corresponding, in the case of a flat-spectrum source at {z) = 2.75, to a B-band luminosity of about O.IL*, 
or a total star formation rate (Salpeter IMF) of ^ 1 Mq yr~^. The brightest has Veoe = 23.6, and is forming 
stars at a rate of ~ 20 yr~^. 

It is interesting to compare the comoving space density of the HDF F300W dropouts brighter than 
^606 = 25 , 7.3 X 10~^ Mpc~^ at (z) = 2.75, with that derived from Steidel et al. (1996a) ground-based 
statistic, 3.6 x 10~*Mpc~^ to 7?. < 25 at redshift {z) = 3.25. Within the errors, the two estimates appear in 
good agreement with one another, especially after accounting for the fact that one probes ~ 0.3 mag fainter 
in the galaxy luminosity function at the HDF average redshift. 

4-2. F45OW Dropouts: Limits to the Luminosity Density at 3.5 < z < 4.5 

In a similar manner, we have identifled 14 F450W dropouts which satisfy the criteria established for 
star-forming galaxies at 3.5 < z < 4.5. Their comoving number density at (z) = 4.0 is 1.5 x 10~'^Mpc~^ 
(2.4 X 10~^ Mpc~'^ for qq = 0.05). From the observed I814 magnitudes, we derive a specific (comoving) 
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emissivity at 1630 A of 

P1630 ~ 5.0 X 10^^ ergs s"^ Hz"^ Mpc"^ (9) 

(2.8 X 10^^ ergss"^ Hz'^ Mpc'^ for go = 0.05). The faintest F450W dropout has a magnitude of /814 = 27.5, 
corresponding, in the case of a flat-spectrum source at {z) = 4, again to a B-band luminosity of about 
O.IL*, while the brightest has J814 = 25. Our samples of C/300 and B450 dropouts therefore reach comparable 
depths. We stress that, contrary to the case at 2 < z < 3.5, we have yet no spectroscopic confirmation of 
the efficiency of our color- selection criteria at 3.5 < z < 4.5. 

The source catalogs of the f/aoo and B^^q dropouts are presented in Tables 2 and 3. For each galaxy we 
report the FOCAS catalog entry number from Williams et al. (1996) (ID), the right ascension and declination 
{oi,5), and the AB isophotal magnitudes in the four HDF passbands.^ The UV dropouts 2-449.0, 3-550.0, 
4-555.1, and 4-676.0 have been spectroscopically identifled by Steidel et al. (1996b) as galaxies at z =2.845, 
2.775, 2.803, and 2.591, respectively. 

5. COSMIC STAR FORMATION HISTORY 

The results of the Canada-France Redshift Survey (CFRS, Lilly et al. 1995) over the redshift range 
< z < 1, the spectroscopic confirmation by Steidel et al. (1996a,b) of the existence of a substantial 
population of luminous, star-forming galaxies at 2.5 < z < 3.5, together with our statistical constraints on 
the redshift distribution of faint galaxies in the Hubble Deep Field, as derived above for 2 < z < 4.5, have 
considerable implications for our understanding of the global history of star and structure formation in the 
universe. Here, we review some of the consequences of these observations, focusing on what can be learned 
about galaxy evolution at early (and late) epochs from integrated quantities over the entire population, 
rather than from a detailed study of individual objects. 

5.1. The Metal Production Density of the Universe 

The UV continuum emission from a galaxy with significant ongoing star formation is totally dominated 

by short-lived massive stars, and is therefore nearly independent of the galaxy history. Moreover, the (rest- 
frame) radiation fiux below 3000 A is a very good measurement of the instantaneous ejection rate of heavy 
elements {Z > 6), since both are directly related to the number of massive stars (Cowie 1988; Songaila, 
Cowie, & Lilly 1990; Madau & Shull 1996): the same stars with m > 10 Mq that manufacture and return 
most of the metals to the ISM also dominate the UV light. The supernova event from a star of 25 M.q injects 
about 4.5 Mq of metals (Woosley & Weaver 1995). At the end of the C-burning phase ~ 17 Mq have been 
converted into helium and carbon (Maeder 1992), with a mass fraction released as radiation of 0.007. For 
each 1 Mq of metals ejected, to first approximation we then expect (17 x 0.007/4.5) = 0.025 M© c^ of energy 
to be radiated away. 

What we are interested in here is the universal rate of ejection of newly synthesized material per unit 
comoving volume. In the approximation of instantaneous recycling, the metal ejection rate (MER) per unit 
volume can be written as (cf Tinsley 1980) 

Pz = '4' J mpzja(l>{m)dm, (10) 

where ip is the star formation rate (SFR) density, 0(m) is the IMF (normalized through the relation 
/ m0(m)dm = 1), and p^m is the stellar yield, i.e., the mass fraction of a star of mass m that is converted 

^ The use of a more conservative criterion which assigns a 2a lower limit to the U300 — -B450 colors of 
galaxies undetected in the F300W bandpass would decrease the number of ultraviolet dropouts to 47, thereby 
decreasing the integrated emissivity given in eq. (8) by a factor of 1.3. The same 2a color criterion applied 
to galaxies undetected in F450W decreases to only 9 the number of blue dropouts, resulting in an emissivity 
which is a factor 1.4 lower than that given in eq. (9). 
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to metals and ejected. The dot denotes differentiation with respect to cosmic time. At short wavelengths, 
the luminosity density radiated per unit frequency during the main sequence phase is related to by 

= O.OOTc^i/j J mfue{'m)fu{m)4){m)dm, (11) 

where fueifn) is the mass fraction of hydrogen burned into helium, and fuinT-) is the normalized spectrum 
of stars of mass m on the main sequence. The detailed conversion from UV luminosity to MER, including 
all evolutionary stages, can be obtained from spectral synthesis and metal-yield calculations. For simplicity, 
we shall ignore here the contributions of planetary nebulae (Renzini &: Voli 1981), Type la supernovae 
(Thielemann, Nomoto, & Yokoi 1986), and stellar winds (Maeder 1992) to the chemical yields, and focus 
on the nucleosynthetic enrichment by supernova explosions from massive stars. We tisc the evolutionary 
models of Bruzual & Chariot (1993) and the Type II stellar yields tabulated by Wooslcy & Weaver (1995) 
(see also Sutherland &; Shull 1996). For a Salpeter IMF including stars in the 0.1 < M < 125 M© mass 
range, a constant SFR, a galaxy age in the interval 0.1-1 Gyr, and solar metallicity, we derive 

P1500 = 4.4 ± 0.2 X 10^9^2 ergss"^ Rz'^ Mpc'^ (12) 

at 1500 A, and 

/92800 = 2.8 ± 0.3 X 10'^^ pz ergs s"^ Hz"^ Mpc"^ (13) 

at 2800 A, where pz is the MER in Mq yr~^ Mpc~ \ The continuum spectrum is fairly flat and is dominated 
by early-type stars. We shall use the relations given above in the subsequent discussion, while noting that the 
conversion efficiency is fairly insensitive to the assumed IMF (as long as the stellar population extends as a 
power-law to massive stars, 50-100 Mq), since the increased metal yield from high mass stars is compensated 
for by a similar increase in the production of UV photons. The advantage then of deriving from the observed 
UV luminosity density a rate of metal ejection, rather than a SFR, is that the latter is instead a sensitive 
function of the IMF slope. To be quantitative, while for a Salpeter IMF the mean yield of returned metals 
is 

j mp^rn(l){rn)dm = 2.4%, (14) 

for a Scalo (1986) IMF - less rich in massive stars - in the same mass range this factor is about 3.3 times 
lower. Hence large errors in the relative rates of star formation at different epochs may occur if the IMF 
varies with cosmic time. Because of these considerations, unless otherwise stated we shall only quote metal 
ejection rates in the following. Note, however, that the real uncertainties in the UV-to-metal conversion 
factors due to age, IMF, metallicity, and population synthesis model well exceed the errors quoted above. 

To relate the amount of heavy elements observed today in various stellar populations and in the gas 
phase to the UV luminosity of distant, star-forming galaxies, we may use the conservation of metals, which 
implies that the sum of the heavy elements stored in stars, Z^p^, and in the gas, Zgpg, is equal to the mass 
of metals ever ejected, i.e., to the integral of equation (10) over cosmic time. An estimate of the baryonic 
mass in galaxies can be obtained by multiplying the observed luminosity density of the local universe by a 
mass-to-light ratio. Since the gaseous content of nearby galaxies is observed to be negligible compared with 
the stellar content (e.g., Rao &; Briggs 1993), the cosmological mass density of heavy elements at the present 
epoch is given by 

M 

pz{0) « Z,p,(0) = Z,pb(0)(— ), (15) 

where Pb(0) is the local blue light density of field galaxies, and M/Lb is the mass-to-blue light ratio of 
visible matter. There are several recent determinations of the B-band luminosity density from large redshift 
surveys. The values obtained for log pb{0) in units of L© Mpc~'^ are 7.98 (Efstathiou, Ellis, & Peterson 
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1988), 7.83 (Loveday et al. 1992), 8.0 (Marzke, Huchra, & Geller 1994), 7.93 (da Costa et al. 1994), and 8.0 
(Ellis et al. 1996). A simple average of all these determinations gives the adopted value: 

/9b(0) = 9.0±1.4x lO^LoMpc-^ (16) 

The mass-to-light ratios for the visible parts of galaxies can be estimated from stellar population modeling. 
The results (in solar units) are in the range 2 < M/Lb < 6, the smaller value referring to a typical 
Population I system like the solar neighborhood (e.g., Kuijken &; Gilmore 1989), the larger value to the 
old stellar populations of elliptical galaxies (e.g., van dcr Marel 1991). Since about 30% of the total blue 
luminosity density is in ellipticals (Efstathiou et al. 1988), we may adopt the luminosity-weighted mean 
value of (M/ Lb) = 3 to derive 

/9z(0) « 5.4±0.8 X 10^(^)MoMpc-=^ (17) 

(of Cowie 1988). While representative stellar metallicities in nearby himinous spirals and ellipticals are 
approximately solar. ~ Zq = 0.02, irregular, blue compact, and dwarf spheroidal galaxies are known to 
be metal poor, > O.IZq (e.g., Pagel & Edmunds 1981). 

The real uncertainty on the metal mass density is difficult to estimate, but is unlikely to exceed a 
factor of 1.5. Although a baryonic mass several times larger than the luminous mass may be present in the 
Galactic halo, metal-rich halo material would be mixed into and over-enrich the disk. Hence, if a substantial 
amount of metals are missing from our census, they are most likely hidden in low-surface brightness galaxies 
(Ferguson & McGaugh 1995), or intergalactic and intracluster gas. 

It is useful at this stage to define a fiducial metal ejection density, 

Sz = Pz{0)tH{0)-^ ~ 4.2 X 10"^ M0 yr-i Mpc'^, (18) 

given by the present-day mass density in heavy elements divided by the current age of the universe. If pz 
observed at a given redshift z is much less than S^, and a large fraction of the luminous baryons observed 
today were already locked into galaxies at this epoch, then either galaxies have already exhausted their 
reservoirs of cold gas or there must be a mechanism which prevent the gas within virialized dark matter 
halos to radiatively cool and turn into stars. By contrast, a metal ejection rate much greater than the fiducial 
value implies that the conversion of gas into stars can be extremely efficient at times. Prom a comparison 
between the observed pz and Sz we may then be able to distinguish between a more or less uniform (with 
cosmic time) SPR and a discontinues one, consisting of a series of bursts of suitable duty cycle. 

According to Gallego et al. (1995), the Ha luminosity density of the local universe is pHa(O) = 1.3 ± 
0.6 X 10'^^ ergss"^ Mpc~^. This value implies^ a MER at the present epoch of 

pz{0) f« 1.1 ± 0.5 X 10"^ M© yr"^ Mpc'^^. (19) 

Equations (18) and (19) suggest that the present rate of production of heavy elements is too low to yield 
the observed element abundances in the Hubble time, and that either star-forming galaxies were much more 
numerous in the past, or individual galaxies must, on average, have passed through a significantly brighter 
phase. 

^ In case-B recombination theory, the Ha luminosity density can be related to the emission rate of ionizing 
photons per unit cosmic volume, niom according to pHa ~ 0.45/iJ^HQ'^ion (Osterbrock 1989). Population 
synthesis galaxy spectra (Bruzual & Chariot 1993) with constant SFR yield the following approximate 
relation between pisoo and nion: ?^ion ~ 0.14pi5oo/^, nearly independent of age in the interval 0.1-1 Gyr. 
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From the CFRS faint galaxy survey (Lilly et al. 1995), we now know that the MER per unit mass was 
in fact considerably larger in the past. The comoving metal ejection density at < z < 1 corresponding to 
the observed specific luminosity density at 2800 A (the "LF-estimated" value of Lilly et al. 1996) is 

\ 3.9 

pz(-z) ~ 1-5 ±0.6 X 10-3 f—tfj M0yr-lMpc-^ (20) 

about 1-5 times higher at z ~ 1 than the current value. This strong evolution is associated with galaxies 
bluer than typical Sbc's. The luminosity function of redder galaxies shows little change back to z ~ 1. 

The sample of Lyman-break star-forming galaxies identified by Steidel et al. (1996a) has a comoving 
density of ^ 3.6 x IQ-^Mpc-^ at (z) = 3.25, and an average continuum specific luminosity at 1500 A of 
10 ergs s" -1 A-i. This gives a comoving volume emissivity of 2.7 x 10^^ ergs s ^ Hz ^ Mpc ^, hence a MER 
density 

pz(3.25) « 6.2 X 10"^ Mq yr"^ Mpc-^ (21) 

comparable with the local value but more than 20 times lower than the z = 1 one. 
Finally, our analysis of the HDF images yields, from equations (8) and (9): 

pz(2.75) « 3.6 X 10-"^ Mq yr'^ Mpc-^ (22) 

and 

/9z(4) « 1.1 X lO"'' Mq yr"^ Mpc'^ (23) 

Note that the MER density at (z) = 2.75 is more than a factor of 5 larger than the Steidel et al. (1996a) 
ground-based value at (z) = 3.25, largely because the HDF deep images probe about 1.7 mag fainter into the 
rest-frame luminosity function at these redshifts. Still, the derived luminosity density is significantly smaller 
than the rate at z ~ 1. It is conceivable that our list of candidates in the HDF may be sampling a significant 
fraction of the galaxy luminosity function at these redshifts, and that the contribution to the luminosity 
density from galaxies below our magnitude threshold is small. Strictly speaking, however, the metal ejection 
rates derived above from the HDF data (and the emissivities given in eqs. [8] and [9]) should be interpreted 
as lower limits to the real values, as they only include the fraction arising in the most actively star- forming, 
young and nearly dust-frcc objects. The same considerations also apply to the sample identified by Steidel 
et al. (1996a). We point out that large numbers of faint galaxies at high redshifts are quite plausible in cold 
dark matter-like cosmologies, because of the high space density of low mass dark matter halos predicted in 
any hierarchical theory of galaxy formation (e.g.. White &; Prenk 1991). 

The derived star and element formation history of the universe is depicted in Figure 9. Although the 
star formation densities beyond z = 2 are only lower limits, together with the Gallego et al. and Lilly et al. 
values they seem consistent with the existence of a peak in the cosmic metal production rate in the redshift 
range 1 < z < 2, as predicted by models of the chemical evolution of the damped Lya absorption systems 
(Pei & Fall 1995). 

5.2. Clues to Galaxy Formation and Evolution 

We may at this stage try to establish a cosmic timetable for the production of heavy elements in 
relatively bright galaxies, keeping in mind the inherent uncertainties associated with the estimates given 
above. In an Einstein-de Sitter universe, the total mass density processed into metals over the redshift range 
< z < Zc ~ 1) is, integrating equation (20) over cosmic time, 

Apz{0) = (8.2 ± 2.6 X 10^ Mq Mpc-3)[(1 + z^f-* - 1] = 3.5 ± 1.1 x 10^ Mq Mpc'^ (24) 



Note that the star formation density at z = quoted by Gallego et al. (1995) is derived assuming a 
Scalo IMF, while Steidel et al. (1996a) adopt a Salpeter IMF in their paper. 
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(Contrary to the measured number densities of objects and rates of star formation, the metal mass density 
does not depend on the assumed cosmological model.) If we define two characteristic epochs of star and 
element formation in galaxies, and zz, as the redshifts by which half of the current stellar and metal 
content of galaxies was formed, then a straightforward comparison between equations (17) and (24), together 
with the fact that most of the stars in the inner luminous parts of galaxies are metal rich, imply < ~ 1, 
or in other words that a significant fraction of the current metal content of galaxies was formed relatively 
late, on a timescale of about 8 Gyr. This is comparable with the decay time of star formation, from 5 to 9 
Gyr, required to reproduce the present-day colors of late-type spirals (e.g., Bruzual h Chariot 1993), and 
suggests the possibility that we may be observing in the rcdshift range z = — 1 the conversion into stars of 
gaseous galactic disks. Pure H i disks may be assembled at some higher redshift, and disk gas continuosly 
replenished as a result of ongoing infall from the surrounding hot halo. The case for a "disk epoch" at z ~ 1 
is strengthened by the observations that galaxies of different luminosities and morphological type evolve 
very differently. A common conclusion of several recent deep redshift surveys (e.g., Glazebrook et al. 1995; 
Lilly et al. 1995; Ellis et al. 1996) is that late-type, gas-rich systems show substantial evolution in number 
and/or luminosity from the present to z ^ 1, and that a large fraction of the observed evolution is associated 
with an increase in the mean surface brightness of luminous disks (Schadc ct al. 1995). By contrast, the 
ellipticals and early-type spirals have been remarkably quiescent over the same redshift interval. 

From stellar population studies we also know that about half of the present-day stars - hence metals - 
are contained into spheroidal systems, i.e., elliptical galaxies and spiral galaxy bulges, and that these formed 
early and rapidly. The major arguments come from studies of Galactic bulge color-magnitude diagrams 
(Ortolani et al. 1995) which indicate ages as old as globular clusters, and from studies of the chemical 
abundances in cluster ellipticals and cluster X-ray gas, whose high a-element ratios suggest that both the 
stars and the intracluster gas must have been enriched on timescales of order 1 Gyr or less (Renzini et al. 
1993). If these arguments arc correct (and they are by no means undisputed), then bulges and ellipticals 
must have experienced a bright starburst phase at high z. Where are these protospheroids? 

There are a few circumstantial pieces of evidence supporting the interpretation that the Lyman-break 
galaxies identified at z > 2.5 from ground-based images may represent the progenitors of present-day lu- 
minous galaxies while forming their spheroidal stellar component. From the strength of their interstellar 
UV absorption lines, and within the assumption that their velocity field is dominated by gravity, Steidel et 
al. (1996a) have shown that galaxies at 2; ~ 3 have masses comparable to that of present-day galaxies. 
Smaller masses would result if interstellar shocks, local to the star-forming regions, contribTitc significantly 
to the velocity field. While the HST images of these galaxies show the presence of compact (half-light radii of 
0.2-0.3 arcsec) cores (Giavalisco et al. 1996), from the deeper HDF data is seems that the F300W dropouts 
may be better described as asymmetric or irregular objects (Abraham et al. 1996; van den Bergh et al. 
1996). 

The observations described above may confirm the rationale that galaxies form stars from the inside 
out, i.e., that spheroidal systems have assembled rather early, far beyond the disk epoch at z ~ 1. At the 
star formation density levels inferred from the HDF images at (z) = 2.75 (see eq. [22]) about 15% of the 
observed mass density of metals at z = would have been formed during the "spheroid epoch" at z > 2. 
This fraction could increase if the star formation is shrouded in dust or takes place in galaxies small enough 
individually to fall below the HDF source detection limit. On the other hand, since the metals we observe 
being formed are a substantial fraction of the entire metal content of galaxies, it appears that star formation 
regions remain largely unobscured by dust throughout much of galaxy formation. 

With the data presented here, two arguments which have implications for models of galaxy formation 
can be made more quantitative. 

• From Efstathiou et al. (1988), the space density of bright ellipticals today is n{> L^) = 2.4 x 
10~^ Mpc~'\ If a significant fraction of their stellar population formed in a single burst of duration 1 Gyr 
early in the history of the universe, a comparable number density of objects should be observed at high-z 
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while forming stars at rates in excess of 50 — 100 Mq yr~^. We find that none of the HDF F300W dropouts 
spectroscopically confirmed by Steidel et al. (1996b) has a SFR in excess of 20 yr"-*^ (Salpeter IMF). 
Assuming a probed redshift interval between 2 and 3.5, we can set an upper limit to the comoving density 
of very high star-forming galaxies at {z) = 2.75, n(> 2OM0yr~^) < 6 x 10~^Mpc~^. Hence there appears 
to be a deficit of very bright objects relative to the expectations of the standard early-and-rapidly-forming 
picture for spheroidal systems. 

• In hierarchical models of galaxy formation, small objects form first and merge together to make larger 
ones. While much of the activity associated with star formation and galaxy merging occurs at relatively 
low redshifts (White & Frenk 1991), there is no period when bulges and ellipticals form rapidly as single 
units and are very bright: rather, galaxies become progressively less luminous, more numerous, and more 
compact at earlier epochs. A common prediction of hierarchical models appears to be the steepening with 
lookback time of the luminosity function due to the large increase of the abundance of objects forming stars 
at modest rates (e.g. Cole et al. 1994). It is then interesting to compare, as a function of redshift, the ratio of 
the space densities of galaxies forming stars at rates in excess of, say, 1 and 10 Mq yr~-^. From the Gallego 
et al. (1995) Ha survey, we estimate n{> lM0yr~^)/n(> lOM0yr~^) k, 50 at the present epoch. A similar 
value is found to characterize, at {z) = 2.75, the distribution of SFRs of the F300W dropouts in the HDF. 
Thus, unless a substantial fraction of the modest star-forming galaxies at these early times are red either 
because are relatively old or because are reddened by dust, there seems to be little evidence in our sample 
for a large enhancement of their space density relative to the abundance of high star-forming objects. 

5. 3. Chemical Enrichment of the Intergalactic Medium 

There is one piece of evidence pointing towards a metal ejection density at z > 3 which is at least as 
large as derived in equation (22) , namely the observations of a significant cosmological mass density of metals 
associated with QSO absorption systems at high-z. We shall focus here on the numerous Lya forest clouds. 
Recent spectra at high S/N and resolution obtained with the Keck telescope (Tytler et al. 1995; Cowie et 
al. 1995) have shown that 50%-60% of the Lya clouds with logA'^ni > 14.5 have undergone some chemical 
enrichment, as evidenced by weak, but measurable C iv absorption lines. The typical inferred metallicities, 
■^iGMi range from 0.003 to 0.01 of solar values, subject to uncertainties of photoionization models. Because 
the Lya clouds have large filling factors and contain a significant fraction of the baryons in the universe, the 
new data are strong evidence for a widespread distribution of metals in the IGM. These metals may have 
been produced in situ or in dense gaseous regions of galaxies; the metal-enriched gas was then expelled from 
the regions of star formation to large distances consistent with the sizes of the Lya clouds. 

Associated with the latter type of enrichment a characteristic MER density can be estimated as 

...mO, . (4 . 10- Me ^- MPO-, (^) (f ) . (.5, 

where JIigm is the baryonic density parameter of the Lya cloud phase, /inj is the fraction of heavy elements 
injected into the IGM during a timescale At, and in is the Hubble time at z = 3. If QSO absorption systems 
trace the bulk of star formation occurring in galaxies at high redshifts (Lanzetta, Wolfe, Sz Turnshek 1995; 
Pei & Fall 1995; Madau & Shull 1996), it is therefore plausible that a population of dwarf galaxies, actively 
forming stars at z > 3 at a rate comparable to the limit inferred from the HDF data, may be responsible 
for the contamination of the IGM at high redshifts. Such galaxies have shallow gravitational potential wells, 
thus allowing easy ejection of the heavy elements produced during their active star-forming phase. 

5.4- The Ionizing Background at z > 4 

In the last few years growing evidence has accumulated on the existence of a significant decline in the 

space density of bright QSOs beyond z ^ 3 (Schneider, Schmidt, & Gunn 1994; Warren, Hewett, &; Osmer 
1994; Kennefick, Djorgovski, & de Carvalho 1995). This shortfall, estimated to be about a factor of 5, 
poses serious problems to the idea that a quasar-dominated UV background is responsible for maintaining 
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the intergalactic gas in a highly ionized state above z ~ 4 (e.g., Haardt & Madau 1996). It may be of 
some interest, in this context, to speculate on the possible contribution of our population of star-forming 
"spheroids" to the metagalactic flux at early epochs (see also Songaila et al. 1990). A minimum value for 
the intensity of the ionizing background at z = 4 is given by 

Jl{4) = ^nH,crit(4)J^iGM ~ 2.4 x lO'^^ ergscm'^ 8"^ Hz"! sr'S (26) 

where nn.crit is the closure hydrogen density. This equation reflects the underlying physics of the ionization 
process: ignoring recombinations, at least one UV photon is required per hydrogen atom to photoionize 
the IGM. The ionizing radiation flux is largely local, as sources at higher redshifts are severely absorbed 
by intervening clouds. At z = 4, the attenuation length is only Al ~ 13Mpc (Madau 1992). The required 
(proper) volume emissivity at 1 ryd is then AttJl/AI. Collecting, and converting as usual an ionizing 
luminosity density into a MER per unit comoving volume, we obtain 

Pz,ion(4) « (2 X 10-^ Me yr"^ Mpc'^) ' , (27) 

where now /esc is the escape fraction into the IGM of Lyman-continuum photons. A comparison between 
equations (23) and (27) provides some rationale for the hypothesis that a population of galaxies at z ~ 4, 
producing metals at a rate which is a few times higher than the lower limit derived from the HDF data, 
might rival quasars as a source of photoionization of the IGM. Larger star formation rates would be required 
if a significant fraction of the UV radiation emitted from stars cannot escape into the intergalactic space, 
/esc < 1- 
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FIGURE CAPTIONS 

Figure 1: (a) Mean cosmic transmission for a source at = 2.5,3.5, and 4.5 {solid lines), as a function 
of observed wavelength. The characteristic staircase profile is due to continuum blanketing from the Lyman 
series. Also plotted are the response functions of the four broad passbands, F300W {dotted line), F450W 
{short-dashed line), F606W {long-dashed line), and F814W {dash-dotted line) used for the Hubble Deep 

Field, {b) Magnitude increments AU300 {dotted line), AB450 {short-dashed lines), AVooe {long-dashed line), 
and A/814 {dash-dotted line), derived by integrating the mean cosmic transmission over the corresponding 
bandpass, as a function of the emission redshift. 

Figure 2: (a) Apparent synthetic colors of a star-forming galaxy with constant SFR (age 0.3 Gyr) plotted 
as a function of the emission redshift. Each pair of curves depicts, from top to bottom, the attenuated colors 
and the colors observed in the case of negligible intergalactic absorption. Solid lines: — -B450. Dashed 
lines: -B450 — Veoe- Dotted lines: Vgoe —-^814- {b) The colors of the four brightest starburst galaxies observed 
by HUT (see text for details). Empty dots: U300 — -8450- Filled dots: -B450 — Veoe- Triangles: Veoe — -f8i4- 
All spectra have been reddened by intergalactic absorption. 

Figure 3: [/300 ~ -B450 vs. ^450 — Vqqq for model galaxies. A total of 103879 synthetic spectra of galaxies 
representing a wide range of ages, star formation histories, metallicities, dust contents, and redshifts, were 
folded through the HDF bandpasses. The galaxies shown as small points are at redshifts less than 2 or 
redshifts greater than 3.5. Galaxies shown as large solid circles arc those in the redshift range 2 < z < 3.5 
with ages less than 10^ yr and extinctions As < 1. Large x's are galaxies in the same redshift range that have 
ages greater than 10^ yr or > 1. Our color selection criteria are shown as the large polygon. Galaxies 
within the polygon in the HDF observations are selected as likely candidates for 2 < 2; < 3.5 objects. The 
selection criteria are U300 — B450 > 1.3, U300 — -B450 > -B450 — -^814 + 1-2, and B450 — -f8i4 < 1-5. 

Figure 4: Selection efficiency for F300W dropouts. This figure provides a measure of color selection in 
isolating galaxies with 2 < 2; < 3.5. The solid histogram shows the fraction of models in each redshift interval 
{Az = 0.2) that meet the color selection criteria. The dotted histogram shows the fraction of models with 
ages less than 10^ yr and extinctions Ab < 2 that meet the selection criteria. The color selection criteria 
are extremely efficient at identifying relatively unobscured star-forming galaxies with 2 < z < 3.5, and 
separating them from the chaff at other redshifts. 

Figure 5: -B450 — Vqqq vs. Vgoe — -^814 for model galaxies. The galaxies shown as small points are at redshifts 
less than 3.5 or redshifts greater than 4.5. Galaxies shown as large solid circles are those in the redshift 
range 3.5 < z < 4.5 with ages less than 10^ yr and extinctions Ab < 1- Large x's are galaxies in the same 
redshift range that have ages greater than 10^ yr or Ab > 1- Our color selection criteria are shown as the 
large polygon. Galaxies within the polygon in the HDF observations are selected as likely candidates for 
3.5 < 2; < 4.5 objects. The selection criteria are B450 — Veoe > 1-7(V606 — -^814) + 0-7, S450 — Veoe > 1-5, 
-B450 - Veoe < 3.5(^606 - -^814) + 1-5, and Veoe - 4i4 < 1-5. 

Figure 6: Selection efficiency for F450W dropouts. This figure provides a measure of color selection in 
isolating galaxies with 3.5 < z < 4.5. The solid histogram shows the fraction of models in each redshift 
interval (A2; = 0.2) that meet the color selection criteria. The dotted histogram shows the fraction of models 
with ages less than 10^ yr and extinctions Ab < 2 that meet the selection criteria. The primary soTirce of 
contamination for the sample is galaxies at slightly lower redshifts with intrinsically red colors. The dashed 
line shows the fraction of the models with ages greater than 10^ yr or > 2 that meet the selection critera. 
If such red objects at 2 < 2; < 3.5 are very common, then our derived metal formation rate at higher redshifts 
may be an overestimate. However, we suspect that such old or dusty interlopers are much less common than 
blue star- forming galaxies at slightly higher redshift. 

Figure 7: C/300 ~ -B450 vs. -B450 — ^814 color-color plot of galaxies in the Hubble Deep Field satistying 
magnitude ranges given in the text, (a) Objects undetected in F300W (with signal-to-noise <1) are plotted 
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as triangles at the la lower limits to their J/soo ~ -6450 colors. Symbol size scales with the I^m magnitude 
of the object. The dashed lines outline the selection region within which we identify candidate 2 < z < 3.5 
objects. (6) Location in the color-color diagram of the more than 60 galaxies in the HDF which have known 
Keck/LRIS spectroscopic redshifts from Steidel et al. (1996b), Cohen et al. (1996), and Moustakas et al. 
(1996). 

Figure 8: S450 — Vqqq vs. Vgoe ~ -^814 color-color plot of galaxies in the Hubble Deep Field satisfying 
magnitude ranges given in the text, (a) The symbol types and sizes are as described for Fig. 7a. The 
selection area we define for candidate 3.5 < 2; < 4.5 galaxies is outlined with dashed lines. (6) Same as Fig. 
7b. 

Figure 9: Element and star formation history of the universe. The data points from various surveys provide 
a measurement or a lower limit to the universal metal ejection density, pz, as a function of rcdshift. For 
a Salpeter IMF, to translate pz into a total star formation density, p*, a factor of 42 should be applied. 
Triangle: Gallego et al. (1995). Filled dots: Lilly et al. (1996). Diagonal cross: lower limit from Steidel 
et al. (1996a). Filled squares: lower limits from the Hubble Deep Field images. The dashed line depicts 
the fiducial rate, E^, given by the mass density of metals observed today divided by the present age of the 
universe (see text for details). A flat cosmology with qq = 0.5 and Hq = 50kms~^ Mpc~^ has been assumed. 
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27.57 


27.47 
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25.98 
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12:36:52.697 
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26.80 


3-231.1 
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>29.20 


>29.95 


27.53 


26.60 


3-512.0 


12:36:56.031 
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>28.50 


27.58 


25.91 
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3-534.2 
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3-675.0 
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>28.67 
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62:12:16.12 
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>29.19 
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4-280.0 
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62:12:31.15 
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29.25 


27.17 


26.88 


4-439.0 


12:36:43.753 


62:12:42.63 


>28.63 


>29.39 


26.06 


24.97 


4-702.0 


12:36:43.540 


62:11:51.82 


28.51 


>29.81 


27.53 


26.65 



2 



1 



0.8 



o 

0.6 

ui 
o 

s 



o 

o 
CO 
u 



0.4 



0.2 











1 




HDF WF1-3 



1 1 1 1 1 1 1 1 1 1 1 1 — i — I 1 1 1 1 1 r 



2.591 



2.803 



2.845 3.226 
,775-lg.§g2 



0.32 



0.503 
0.562 



0.96 



0.299 



0.06913 



?^9 



-1 







1 

^450 ~ ^814 



HDF WF1-3 



4 - 



1 1 1 r 



3 - 



to 
o 

CO 



o 
QQ 



- 



-1 



-1 



2. 



1 1 1 — 7~i 1 1 r 

/ 



/ 



1 1 1 1 1 1 r 



/ 



/ 



/ 



0.6i 



^ _ _ _ / °-S!^2 0, 



3.226 0.96 
2.803 



0.37 



32 




773 .355 



0.958 







1 

V - I 

^606 ^814 



-2.5 



1 — r 



1 — 1 — r 



1 — r 



1 



-3 - 



CO 

O 

a 



>;-3.5 



A 



1.5 o 



2 



•-J 



•Q. 

o 



A 



A 

■4 



-2.5 



o 

I 

u 



-4.5 



J I I L 



J I I L 



1 



J I I L 



1 



J I I L 



J I I L 







2 3 
redshift 



